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Polymer ejection from strong spherical confinement 
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We examine the ejection of an initially strongly confined flexible polymer from a spherical capsid through 
a nanoscale pore. We use molecular dynamics for unprecedentedly high initial monomer densities. We show 
that the time for an individual monomer to eject grows exponentially with the number of ejected monomers. 
By measurements of the force at the pore we show this dependence to be a consequence of the excess free 
energy of the polymer due to confinement growing exponentially with the number of monomers initially inside 
the capsid. This growth relates closely to the divergence of mixing energy in the Flory-Huggins theory at large 
concentration. We show that the pressure inside the capsid driving the ejection dominates the process that is 
characterized by the ejection time growing linearly with the lengths of different polymers. Waiting time profiles 
would indicate that the superlinear dependence obtained for polymers amenable to computer simulations results 
from a finite-size effect due to the final retraction of polymers’ tails from capsids. 

PACS numbers: 87.15.A-,82.35.Lr,82.37.-j 


I. INTRODUCTION 

Processes involving macromolecules in confinements are 
studied intensely due to their significance in biology and po¬ 
tential technological and medical applications. An important 
class of such processes is the capsid ejection where a poly¬ 
mer is initially in a compact conformation inside a capsid and 
then ejects outside through a pore a few nanometers wide. 
By far the most important biological process of this class is 
the viral packaging in and ejection from bacteriophages [1- 
9]. The ejection of double-stranded (ds) DNA is clearly the 
most studied case. Ejection processes of biopolymers such as 
RNA and single-stranded (ss) DNA, although less studied, are 
found in vivo and are biologically highly relevant. The recent 
technological advancement in using engineered viral capsids 
for drug delivery [10] further emphasizes the need to under¬ 
stand the fundamentals of biopolymer packaging in and ejec¬ 
tion from a capsid. While dsDNA is modeled as a semiflexible 
polymer, ssDNA and RNA among many other biopolymers 
are fully flexible [11]. 

There is a wide gap between the conditions under which 
biopolymers eject from viral capsids found in vivo and those 
taken to prevail in computer simulations and theoretical anal¬ 
ysis. In the bacteriophages the viral DNAs and RNAs are 
packed to almost crystalline densities [3], realization of which 
is beyond most computer simulations. Computational investi¬ 
gations of some specific characteristics using close-to-realistic 
model polymers packed to high densities are typically done 
using some form of probabilistic Metropolis sampling, see 
e.g. [12]. On the other hand investigations using dynamically 
more realistic molecular dynamics (MD) based methods typ¬ 
ically aim at detailed modeling of a specific polymer and are 
not very conclusive with regard to general characterization of 
the polymer ejection. Due to restrictions placed by the de¬ 
tailed model on polymer lengths and on statistics no attempt 
is made to relate to the existing theoretical formulations. 

Our motivation is to determine how applicable the avail¬ 
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able theoretical formulations to the real-world ejection pro¬ 
cesses are. The available theories do not explicitly take into 
account the bending rigidity in the case of dsDNA. Still, they 
are commonly used as a basis for understanding also the ejec¬ 
tion of dsDNA. The current understanding is that a packed 
dsDNA assumes a spool-like conformation inside the capsid. 
However, in the theoretical formulations polymers start from 
disordered conformations. Blob-scaling picture that describes 
polymers under semi-dilute conditions [13] was used in the 
derivations. This, in effect, means that the theories are not 
guaranteed to be valid for very strongly confined polymers 
found in vivo. 

Here, we set out to study the ejection of fully flexible poly¬ 
mers from spherical confinements. The studied confinements 
are clearly stronger than those in previous studies but still 
much weaker than those found in vivo. The flexible freely- 
jointed chains (FJC) start from random initial conformation, 
so that close comparison with the theoretical predictions can 
be made. In short, we use a generic polymer model but realis¬ 
tic MD based dynamics to characterize the ejection dynamics 
of fully flexible polymers initially under very strong spherical 
confinement. The effect of polymer bending rigidity and the 
accompanying spooled initial conformation that are important 
in the case of dsDNA will be reported in a future publication. 

Due to the computationally effective implementation of the 
capsid geometry sufficient statistics can be obtained for poly¬ 
mer ejection processes starting from the high initial monomer 
densities. In our capsid model boundary conditions are im¬ 
posed instead of introducing a repulsive potential between 
polymer and the capsid. This way external forces that could 
affect the ejection dynamics are avoided. Accordingly, we can 
obtain the polymer’s excess energy due to confinement. We 
show that the form of this energy, which decidedly deviates 
from the predictions based on the blob-scaling picture widely 
used for characterization of confined polymers, determines the 
ejection dynamics. 

We model Brownian heath bath with stochastic rotation dy¬ 
namics (SRD) where also hydrodynamics can be included. In 
order to keep the analysis simple, we use the method here 
with hydrodynamic al interactions switched off. The effects 
of hydrodynamics will be analyzed in detail in the forthcom- 
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ing paper. In what follows, we first outline the understand¬ 
ing previously obtained from the blob-scaling picture. We 
then describe our simulation method after which the results 
are presented and analyzed. Lastly, we summarize and recap 
the main conclusions. 


II. CURRENT THEORETICAL UNDERSTANDING 


The evolution of the theoretical understanding of the ejec¬ 
tion of a fully flexible polymer from a capsid was initiated by 
an investigation by Muthukumar where Monte Carlo simula¬ 
tions were used [1]. Assuming the excess energy due to con¬ 
finement to be AF/ksT ^ Nq/Rq^" a scaling prediction for 
the ejection time of the form r ~ = Ao(Ao/po)^^^^ 

was obtained, where Nq is the degree of polymerization, Rq 
the capsid radius, po the initial monomer density in the cap¬ 
sid, and ly the Flory exponent. This scaling was confirmed 
by MC simulations. In a later work [7] the excess energy 
due to spherical confinement was taken to follow the scal¬ 
ing AF/keT ~ ~ where 

00 = Noa^ /Rq is the initial monomer volume fraction and a 
is the monomer length. This scaling law, first introduced by 
Grosberg and Khokhlov [14], was shown to be correct for the 
spherically confined polymer at 0o where the blob scaling in 
the semidilute conditions is valid [15]. This led to the scaling 
relation r ~ again confirmed by MC simu¬ 

lations. The different regimes and short ranges explain why 
different scaling relations in these studies were corroborated 
by simulations. 

The unified framework [8] for polymer decompression pro¬ 
cesses presented a fairly complete view of the ejection pro¬ 
cess. Assuming uniform polymer conformation the excess 
confinement energy inside a spherical capsid was derived in 
the same form as in [14], 

AF/ksT ^ , (1) 

where N (t) is the number of monomers inside the capsid. 

The resultant driving force ~ kBT/^{t) was taken to be 
exerted on the monomer residing at the pore. The overall 
dissipation takes place close to the pore, within the range 
of the correlation length ^{t), where there exists a veloc¬ 
ity gradient of segments of the order ^ aN{t)/^{t). Ac¬ 
cordingly, the dissipation term was evaluated as TS{t) = 
r][N(t)/^{t)Y. The excess confinement energy is dis¬ 
sipated at the rate AF{t) = —TS{t), from which the time 
evolution was obtained as 

N{t)=No{l + t/Tif, (2) 

where the exponent /3 = (1 — 3r^)/[2(l — v)] and the 

time constant ri ~ ^Nq. tq ~ r/a^/(ksT) 

is the monomer scale time constant, where rj is the solu¬ 
tion viscosity. For the pressure-driven part the scaling r ^ 
to0q was obtained. [8] 



FIG. 1. (Color online) Snapshots of a simulated polymer ejection. 
(Images created using VMD [21] and POV-Ray [22]) 


Blob-scaling assumption was used in the theories outlined 
above. Also in the strong confinement, as defined in [8, 15], 
the confined polymer chain is assumed reminiscent of a 
semidilute polymer solution and the description is based on 
the blob-scaling assumption. 


III. THE COMPUTATIONAL MODEL 


For the simulation of the ejection dynamics we use a hybrid 
computational method where the time-integration of the poly¬ 
mer is performed by MD implemented by the velocity Verlet 
(vV) algorithm [16, 17]. The polymer is immersed in a sol¬ 
vent that is modeled by SRD [18, 19]. Here, we exclude the 
hydrodynamics for the better understanding of the fundamen¬ 
tal ejection dynamics and more straightforward comparison 
with the existing understanding based on the blob picture. The 
spherical capsid is modeled as a shell with rigid walls impos¬ 
ing slip and no-slip boundary conditions for polymers and sol¬ 
vent, respectively. The pore is modeled as a cylindrical hole in 
the shell. The simulation geometry and snapshots of an eject¬ 
ing polymer are depicted in Fig. 1 . The polymer bond length 
at rest is approximately 1. The pore radius Rp is 0.4 for the 
polymer and 0.8 for the solvent. The thickness of the capsid 
wall is 3. The capsid geometry was created using computa¬ 
tionally effective constructive solid geometry technique [20]. 


The polymer is modeled as a chain of point-like beads with 
mass mb. Adjacent beads are connected via the FENE poten¬ 


tial Up = 1 - 




where r 


is the distance between adjacent beads and K and rmax are po¬ 
tential parameters describing the strength and maximum dis¬ 
tance limit of adjacent beads. A Lennard-Jones potential acts 


between all beads: Up,] = 4.8e 


(*)-(*) 


-F 1.2e 


for Tij < \/2a and Ulj = 0 for > -^a. e and a are po¬ 
tential parameters and is the distance between beads i and 
j. The potential parameters are chosen as cr = 1.0, e = 1.0, 
K = 30/(7^, and r^ax = l-5cr in reduced units [23]. 

We use initial monomer density po = No/{^ttRq) instead 
of volume fraction 0o- Accordingly, the number of monomers 
Nq corresponding to po = 1 is by the factor larger than A"o 
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corresponding to (/)o = 1. On the other hand, (po was used for 
hard spheres, the use of which is not possible in high-density 
MD simulations. A value of po constitutes a slightly lower 
compression in a system using soft potentials than in one us¬ 
ing hard-sphere potentials. This effect is more than compen¬ 
sated by our using po instead of po, so effectively the systems 
simulated here start from more compressed states than those 
dealt with in [1, 7, 8]. In the initial conformations four beads 
are inside the pore so the total length of the polymers is Nq+A. 

SRD solvent consists of point-like particles whose dynam¬ 
ics can be divided into two steps. In the streaming step the 
solvent particle positions are ballistically propagated in dis¬ 
crete time steps. The interactions between particles are taken 
into account in the collision step. Here the random parts of 
the velocities for the polymer and solvent particles divided 
into cubic cells of unit edge lengths are rotated by the angle 
a = 37r/4 around an axis chosen randomly for each cell. In 
the present case of Brownian heat bath velocities are randomly 
exchanged between all particles after the collision step. The 
solvent is kept at the constant temperature of ksT = 1.0 by 
scaling the random parts of particle velocities such that the 
equipartition theorem holds at all times. In order to maintain 
Galilean invariance, the grid is shifted randomly at each time 
step [24]. For our simulations the vV time step St = 0.0002 
and the SRD time step At = 0.5. MD and SRD steps are per¬ 
formed in turns such that after At/St vV steps a single SRD 
step is performed. More detailed accounts of the SRD method 
can be found in [18, 19]. 

The unit of length in simulations is typically mapped to 
real-world length scale in various ways depending e.g. on the 
polymer that is modeled. Here we present a mapping when 
the FJC polymer of our simulation is taken to model a single 
stranded DNA (ssDNA). The persistence length of ssDNA is 
of the order of 4 nm [25]. In our FJC the the persistence length 
is Xp = ^b, where b is the bond length of about 1 in our sim¬ 
ulation units. Accordingly, one simulation unit a corresponds 
to about 8 nm and a polymer of length Nq = 200 to a ss¬ 
DNA approximately 1600 nm long. Since the base in ssDNA 
is about 0.37 nm long [26], a polymer of length Nq = 200 
corresponds to about 4324 bases. In our simulations, the cap¬ 
sid inner radius varies from 1.6 (N = 25, po = 1.5) to 5.1 
{Nq = 283, pq = 0.5), which corresponds to the range from 
12.8 to 40.8 nm. 

To gain understanding about the time scales involved in the 
ejection process, we measured the relaxation times of free 
polymers starting from fully stretched conformations. The 
measured radii of gyration of polymers follow approximately 

Rg{t) ~ [i?g(0) - Rg{oo)] exp (-f/Trelax) + Rg{oo), (3) 

where relaxation times r^eiax = 1642, 8883, and 28780 for 
Nq = 50,100, and 200, respectively. Hence, in our simula¬ 
tions the ejection takes from 4.2 to 12.5 times the relaxation 
time for Nq = 50 and from 2.5 to 7.1 times the relaxation 
time for Nq = 200. In principle, this can be compared to 
relaxation time of a real ssDNA. In Ref. [27] the authors mea¬ 
sured relaxation times for single stranded A-DNA. They find 
that the relaxation time is a = 0.19 s in solution of viscos¬ 
ity 15 cR The relaxation times from experiments and our sim¬ 


ulations are by no means comparable. In experiments polymer 
conformations evolve by Zimm dynamics, of course, whereas 
our simulated polymers perform Rouse dynamics. In addition, 
the length of the A-DNA is 48502 bases long [28], which is al¬ 
most ten time as long as in our simulations. Also the viscosity 
in their experiments is much larger than in our simulations. 


IV. RESULTS 

A. Ejection time and waiting times 

Fig. 2 (a) shows the ejection time r vs number of monomers 
Nq initially in random conformations inside the capsid for 
different initial densities pQ. In accordance with previous 
findings [1-4, 6-9], scaling t ^ Nq is obtained. Here, 
P = 1.362±0.05,1.325±0.04,1.293±0.04,1.299±0.03,and 
1.300 ±0.04 for Po = 0.5, 0.75,1,1.25, and 1.5, respectively. 
(As stated in the Introduction, using initial random conforma¬ 
tions is essential for evaluating the validity of the blob-scaling 
picture. We also checked that for FJC polymers the initial 
conformation does not have a strong influence on the ejection 
dynamics. For the simulated densities po G {0.5,1.0,1.5} 
we obtained essentially the same /3 when the ejections started 
from spooled conformations.) /? decreases slightly with in¬ 
creasing initial monomer density for po < 1 in accordance 
with our previous findings [9] and statement on the case of 
driven translocation when the pore friction increases [29, 30]. 
Most of the effective friction 7 is exerted on the ejecting poly¬ 
mer in the vicinity of the pore through which the pressure 
pushes it. 7 increases with pQ. As 7 increases r ^ Nq, 
P ^ 1. In accordance with this reasoning, for polymers of 
Nq > 50, T oc 1 /Pq is obtained very precisely. 



Nq Rq 


FIG. 2. (Color online) (a) Ejection time r as a function of poly¬ 
mer length No for different initial monomer densities po. The points 
are averages over 40-50 runs. The lines are the fitted curves of the 
form T ~ N^. P = 1.362 ± 0.05, 1.325 ± 0.04, 1.293 ± 0.04, 
1.299 ± 0.03, and 1.300 ± 0.04 for po = 0.5, 0.75, 1, 1.25, and 
1.5, respectively, (b) r as a function of capsid radius Ro for different 
No, T ~ Rq. From bottom to top: = 2.71, 2.77, 2.90, 3.05, 2.97, 
3.23, 3.19, and 2.92 for No = 25, 35, 50, 71, 100, 141, 200, and 
283, respectively. Both figures in log-log scale. 

Fig. 2 (b) shows r as a function of capsid radius Rq. Com¬ 
paring the obtained scaling r ~ Rq, where (/ G [2.71, 3.23], 
with the scaling in Fig. 3 (a) of [ 8 ], where ( decreases from 
5 toward 2 when Ro/a decreases from 10^ toward 1, indi¬ 
cates that the ejection takes place in the very strongly confined 
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regime, that is to say, in stronger confinement than required 
for the strong confinement as defined in [8, 15]. We also mea¬ 
sured the radii of gyration Rg for the polymers’ initial con¬ 
formations inside the capsids. Spherical scaling Rg ^ iVg, 
where 9 k, 0.33, was obtained for all p^, also confirming 
strong confinement. In spite of the very strong confinement 
the ejected polymer segment remains close to equilibrium, 
since we measured Rg{s) ~ ® for the ejected segment at 

different stages. Magnitudes for Rg{s) were also close to the 
equilibrium values. 

Although ejection times r scale with Nq, as seen from the 
endpoints of the cumulative waiting times f(iVo), cumulative 
waiting times t{s) as such do not scale but grow exponentially 
with s, see Figs. 3 (a) and (b). Accordingly, waiting times 
of the ejecting monomers fu,(s) = t{s) — t{s — 1) follow 
the exponential form, ~ exp(Ks), as seen in Fig. 3 (c), 
in accord with the waiting time profile found in [9] for the 
symmetric pore. 



FIG. 3. (Color online) (a) The cumulative waiting time of simula¬ 
tions where po = 1.0. Averages over 50 runs. The times when 
the end monomers of polymers of different lengths escape scale as 
T ~ (solid line). Reaction coordinates at instants when the 

inner part of the polymer exerts no force to the bead at the pore (■): 
sjrl) ~ g]-^y line). Reaction coordinate when num¬ 

ber of beads inside the capsid N = go (k) (b) The cumulative wait¬ 
ing time vs normalized reaction coordinate s/Nq for po = 1.0 and 
different Nq. Squares (■) and triangles (A) correspond to those in 
(a), (c) Waiting times tw and the inverse of the measured force at 
the pore 1// (arbitrary units) for po = 1.0 and po = 1.5 when 
No = 200. V is the capsid volume, (d) Waiting times tw vs. nor¬ 
malized reaction coordinate s/Nq for po = 1.0 and No = 50,100, 
and 200. 


B. Excess energy due to confinement 

The form ~ exp(Ks), or equivalently, tw{N{t)) ~ 
exp{— kN( t)) is at odds with AF based on the blob-scaling 
picture, Eq. (1). To determine how largely AF determines 
tw{s) we measured the force / that has to be exerted on 
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FIG. 4. (Color online) (a) Force at the pore entrance in equilib¬ 
rium as a function of the monomer density inside the capsid po for 
a capsid of volume V — 200/1.50 (No — 200 and po = 1.50). 
The curve that is higher in the range p € (0.3, 1.1) is the best 
fit to the data of the force predicted by the Flory-Huggins the¬ 
ory. The other curve is an exponential curve fitted to the data, (b) 
Number of beads inside the capsid N{t) during the ejection for 
iVo = 25, 35, 50, 71,100, 141, 200, and 283 with fits to Eq. (5) and 
the largest No fitted to Eq. (2). 


the monomer at the pore to keep it in position at different s. 
Fig. 3 (c) shows and 1// for po = 1.0 and 1.5. 1/ f ^ tw 
is seen to hold well for longer chains throughout the ejec¬ 
tion until the start of the final retraction of the remaining seg¬ 
ment in the capsid. For po = 1-5 this relation holds better, 
as is expected due to increased pressure driving the ejection. 
The initial deviation is caused by increased jamming at higher 
densities. Hence, ejection rate is dominantly determined by 
AF{N) at all stages except for the final retraction. 

What is the origin of the exponential dependence / = 
C exp[i3A^(f)]? Generally, it can be stated that it arises due to 
excluded volume interactions as inside the capsid monomers 
interact individually rather than as ensembles of blobs and the 
higher order interaction terms become important. Cacciuto 
and Luijten point out that the blob scaling description breaks 
down at sufficiently high concentration [31]. They proposed a 
different scaling to set in at this concentration due to screening 
effects. 

Indeed, the average monomer number per blob is very low 
at concentrations used in our simulations. The blob radius 
Th = AqNI^, where Ni, is the average number of monomers 
in a blob. Aq = 0.6 was used for the self-avoiding chains 
in [31]. The average density within a blob pt = Nb/{A/ 
has to equal the density inside the capsid, which gives us Nb = 
(4/37rAgp)i/(i-3‘^). Forp = 0.5 Nb^?, and for p = 1 Nb ~ 
1. Since for monomer densities used in our simulations the 
average number of monomers per blob is between 1 and 3, it is 
obvious that the blob-scaling picture cannot be valid for these 
densities, and monomer-monomer interactions determine the 
excess energy due to confinement. 

In SRD the polymer is immersed in a solvent. Accordingly, 
the excess energy due to confinement should be derived start¬ 
ing from the free energy of the mixing of the polymer and 
the solvent. We may approximate this by the mixing free en¬ 
ergy per lattice site for a polymer in Flory-Huggins theory 
= N-^(l)ln{cj)) -F (1 - (^) ln(l - 0) -Fx(l - </'), where 
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X is a constant [32]. Flory-Huggins theory calculates this en¬ 
ergy for a polymer consisting of hard spheres. Here, (p is the 
fraction of lattice sites occupied by polymer segments, fm di¬ 
verges as (1 — (p) ln(l — (p), when cp ^ 1. For the hard-sphere 
lattice polymer the prefactor would go to zero as (p approaches 
1. For our continuum model, where soft potentials are used 
for the polymer the prefactor for the effective ^ oc p is not 
exactly of this form and does not have to vanish for (pef / = 1- 
Fig. 4 (a) shows the divergent term of the Flory-Huggins mix¬ 
ing energy and the exponential form htted to the measured 
force at the pore f.f = C exp[BN (f)] gives the best ht, but 
also the Flory-Huggins form aligns with the data reasonably 
well. We conclude that the mixing entropy explains our mea¬ 
sured free energy due to conhnement. The deviation from the 
exact form derived exactly within the Flory-Huggins theory 
arises due to the soft interaction potentials and the continuous 
space in our model. [33] 


C. Time dependence of the monomer number in the capsid 


Using the measured dependence / = C'exp[^iV(f)/iVo], 
we can derive N{t) in the framework presented in [8]. If the 
monomers are packed inside the capsid by force applied at 
the pore, then the excess energy that results from packing N 
monomers is AF fv fiAk, where fi is the force re¬ 

quired to move the bead i at the inner pore opening into the 
capsid and the bead i -f 1 in its place. Ali is the distance the 
bead i -\- 1 needs to be moved. Since fi is measured for in¬ 
dividual beads, AU = a ~ 1, Vi. In the limit Ali —>■ 0 and 
the minimum necessary force being applied continuously on 
the polymer at the inner pore opening AF = f{n)dn = 

C exp{An/NQ)dn = (CA^oM)[exp(i?iV) — D], where 
C oc ksT, and D are constants. Relating the rate of change 
of this energy to the overall dissipation (see text after Eq. (1)), 
AF{t) = -TS{t) Ri we get 


Nit) 


C AN it)/No 

Va^fit) 


(4) 


For the high monomer densities monomers close to the pore 
are pushed out from the dense initial conformation. Hence for 
the large part of the ejection the correlation length is not ex¬ 
pected to change appreciably. Approximating the correlation 
length to be constant f{t) = f and using the initial condition 
7V(f = 0) = Nq, the solution is given in the form 


N{t) 



AC 

No'qa'^f 



(5) 


Fig. 4 (b) shows the ht of Eq. (5) to the measured Nit). A ht 
of Nit) given by Eq. (2) for Nq = 283 is given for reference. 
Eq. (5) is seen to describe Nit) with excellent precision as it 
should given that t^^ (x 1// almost throughout the ejection. 


D. Ejection time grows linearly with polymer length 

Eor constant pQ, the waiting times as functions of the 
normalized coordinate s/Nq fall on the same curve except for 


the hnal retraction, see Eig. 3 (d), so for this part t ^ Nq. 
The following calculation also conhrms this: Eigs. 3 (a) and 
(b) show the times when the measured force exerted by 
the polymer segment inside the capsid on a monomer at the 
pore is zero. These points, determined by measuring the force 
at the pore needed to keep segments of different lengths N 
completely inside the capsid, obey iV(f = T 2 ) = KNq, where 
K m 0.38. Solving Eq. (5) for f gives 


t; = - e-^) = PNq, P > 0 (6) 

for the pressure-driven part. This linear dependence does not 
result from Nit = T 2 ) ^ Nq but is mainly due to AF not 
scaling with N, which can be shown by replacing AF of 
Eq. (1), by our measured exponential form in the framework 
presented in [8]. This leads to linear dependence for large Nq. 

Erom Eig. 3 (d) it is seen that (s /No ) for different Nq de¬ 
viate from the common form only at the hnal stage of ejection 
when the remaining part of the polymer retracts from the cap¬ 
sid. Retraction speed increases identically for all polymers. 
After T 2 the force at the pore due to entropic imbalance be¬ 
tween the outside and inside of the capsid causes the tension 
to propagate from the pore along the polymer segment in the 
capsid. The tensed segment As grows identically for all Nq. 
Hence, As /Nq will be larger for smaller Nq resulting in re¬ 
traction starting at smaller s/Nq for shorter polymers, as seen 
in Eig. 3 (d). This hnite-size effect results in superlinear scal¬ 
ing T ~ Nq, /? > 1. Accordingly, for asymptotically long 
polymers linear dependence t ^ Nq would be obtained. 


E. Final remarks 

Einally, it is in place to note that to determine the free en¬ 
ergy of the polymer in a very strong spherical conhnement, 
as we have done, it is essential to use the generic EJC poly¬ 
mer model. Else, for example the validity of the blob-scaling 
picture could not be evaluated. EJC is the relevant model for 
example when simulating the ejection of ssDNA or RNA. In 
spite of our generic polymer model the situation is in many re¬ 
spects analogous to a realistic dsDNA in a capsid. The pack¬ 
ing force for the dsDNA in a capsid is mainly determined by 
the long-ranged electrostatic interactions. Hence, in both our 
model and the dsDNA the dominating interactions are not the 
interactions of blobs but interactions of individual monomers 
via repulsive potentials, in our case the Lennard-Jones po¬ 
tential and in the case of a DNA the Coulombic potential. 
Exponential-looking dependencies of the packing force on 
monomer concentration have, indeed, been obtained experi¬ 
mentally, see e.g. [3]. Also in the related computational study 
a dependence resembling exponential was found [34]. How¬ 
ever, since the dependencies were plotted in linear scale and 
the ranges were quite small, dependency in these two studies 
cannot be determined with certainty. A recent computational 
study for semihexible polymers also shows a packaging force 
dependence that does not scale with Nq [35]. 
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V. CONCLUSIONS 

In summary, we have investigated in detail the ejection of 
flexible polymers from spherical capsids through a nanoscale 
pore via computer simulations using realistic dynamics. The 
ejection dynamics and the pertaining excess energy due to 
confinement AF were analyzed via measured waiting time 
profiles and forces exerted on polymers at the pore. We found 
that the waiting times grow exponentially with the number 
of ejected monomers and that the force / measured at the pore 
increases exponentially with the number of monomers in the 
capsid N. We showed that AF then must grow exponentially 
with N, which we addressed to be due to the higher-order 
terms in monomer-monomer interactions. This exponential 
dependence is very similar to the divergence of the mixing 
energy in the Flory-Huggins theory at high concentration. We 
also found that tw ^ 1 // holds well for the simulated den¬ 
sities and that, accordingly, AF determines the ejection dy¬ 
namics for such strongly confined polymers. We showed that 
this AF results in the ejection time r growing linearly with 
the polymer length Nq. The measured superlinear scaling 
T ~ Nq , P > 1, results from a finite-size effect due to the 


final retraction of polymers’ tails from capsids. 

Importantly, for densities that are larger than those asso¬ 
ciated with strong confinement but still moderate compared 
with realistic densities the strong monomer-monomer interac¬ 
tions result in the excess energy due to confinement increasing 
exponentially with the number of monomers in the capsid and 
ejection time increasing linearly with polymer length. This is 
in stark contrast with the previous results on ejection dynam¬ 
ics that are valid only in the semidilute regime where blob¬ 
scaling picture applies. Also of importance is that, contrary to 
some claims, for these densities the ejection easily completes 
without any assisting mechanism such as flow. 
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